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Abstract. We give statistical definitions of the length, I, of a loose prime knot tied 
into a long, fluctuating ring macromolecule. Monte Carlo results for the equilibrium, 
good solvent regime show that (I) ~ N*, where N is the ring length and t ~ 0.75 is 
independent of the knot type. In the collapsed regime below the theta temperature, 
length determinations based on the entropic competition of different knots within the 
same ring show knot derealization (t ~ 1). 
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Long polymers are remarkable systems in which topological entanglement is almost 
unavoidable 13 an d shows consequences pertaining to many different scientific 
contexts. In chemistry the synthesis of macromolecules with nontrivial topology, like 
catenanes, is well established The role of knots and links in the biology of DNA 
and the existence of enzymes like topoisomerases, which control these entanglements, 
have been at the basis of an explosion of interest in topology among biologists, chemists, 
mathematicians, and physicists since the eighties [3J |3] . Most recently, techniques for 
manipulating single molecules opened new possibilities of knot creation and study |EllZj- 

In all these examples the knowledge of the degree of localization of the knots within 
the chains is of primary importance. The action of topoisomerases E] certainly 
depends on how localized the topological entanglement is The folding dynamics of 
a knotted protein [9 J should also be strongly influenced by the tightness of the knot. 
When a knot diffuses along a DNA molecule stretched by molecular tweezers, the length 
of the entangled region is also essential since it determines the diffusion coefficient. If 
the knot is tight enough this length can be approximated by that of the knot in its ideal 
form [7 . An ideal knot is realized with the shortest piece of rope of constant diameter 
and is the only one for which a precise definition of length exists (TQIEJE]- Of course, 
this model is inadequate if the molecule is not well stretched and the knotted region 
fluctuates substantially. 

In spite of several indications that prime knots of polymers in good solvent are 
rather well localized [T3J E], proofs and direct quantitative descriptions of this 
localization are still missing. This is primarily due to difficulties related to the notion 
of a knotted open string [T3]. Indeed, to localize the knot requires to identify an open 
portion of the ring which is "knotted". However, to assert that an open piece of a 
closed rope "contains" a knot, is ambiguous. Mathematically there are no knots in 
an open string since continuous deformations can always bring it into an untangled 
shape. Substantial progress in the description of localization properties of topological 
entanglement has been made for flat knots jTHj. These models describe three-dimensional 
knotted polymers fully adsorbed on a plane with the simplifying feature that the number 
of overlaps is restricted to the minimum compatible with the topology. In this way 
knot length can be defined and studied both numerically and analytically [Uj. It 
was also discovered that flat knots become delocalized when the adsorbed polymer 
undergoes theta collapse [THl EH|. A demonstration of the conjectured jTH] analogous 
derealization phenomenon for compact knotted polymers in three dimensions remains 
an open challenge [20J. 

In this Letter we propose a statistical description of the lenght, I, of a knot within a 
circular fluctuating polymers with N monomers. We find that (I) ~ N* with t ~ 0.75 in 
the good solvent regime (weak localization), and t ~ 1 (derealization) in the collapsed 
phase pi] . 

We consider knotted self-avoiding polygons (SAP) on cubic lattice, which are a good 
model for long polymers in a good solvent [21]. If necessary, a theta collapse (2I| into 
compact configurations can be induced by switching on an attractive interaction between 
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nearest neighbor lattice sites visited not consecutively by the SAP ^H] and by lowering 
enough the temperature T. N is the number of lattice edges occupied by the SAP. A 
BFACF Monte Carlo dynamics [22] preserving the topology of the SAP j2H| enabled the 
sampling of equilibrium configurations of chains of variable N. To increase mobility a 
multiple Markov chains in the space of the edge fugacity was implemented. 

A first strategy of knot length determination is as follows. In each sampled 
configuration, various open portions of the SAP are considered and for each one of them 
a closure is made by joining its ends, A and B, with an off-lattice path (Fig.l). This 
path is chosen in order to minimize the risk of knot modifications or disentanglements 
in the resulting new ring. This risk can not be fully avoided and is a possible source 
of systematic errors. Once the new ring has been constructed, the computation of a 
topological invariant, namely the Alexander polynomial [26J A(z) at z — — 1, allows to 
verify the presence and the type of knot [2Z|- The length / of the knot in a given SAP 
configuration can then be identified with the shortest ring portion still displaying the 
original knot, among a large set obtained by several cutting and closing operations [28J. 
The plot of the average length of a trefoil knot (3i) j2H] as a function of N (Fig. 2) shows 
full consistency with the scaling law 

(I) ~ N* (1) 

with t = 0.74 ±0.14 which is also robust with respect to a change of prime knot type. For 




Figure 1. Once extracted the open string portion (AB) from the SAP, we determine 
its center of mass C. An off- lattice planar 4-edges polygonal is attached to the extrema 
A and B. The polygonal lies on the plane CAB. D,E and F are all chosen far away 
from C in order to avoid additional entanglement with the open string. 
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Figure 2. (I) of the trefoil (3i) knot as a fuction of N. The configurations sampled 
were ~ 10 6 in total. The brackets (.) denote averages taken over data binned around 
fixed values of N. The inset shows the log-log version of the plot. 



example, if we replace the 3i by a 4i or 5i knot j26j . Eq.(l) remains valid with the same 
t within confidence limits. Unfortunately, one can not estimate the possible systematic 
errors arising from the cutting and closing procedure used in this analysis [2E], which 
would be a serious handicap for applications to compact polymers. Thus, one needs to 
validate the law in Eq.(l) by alternative, consistent methods of length determination. 

Consider a SAP partitioned into two loops by a narrow slip link. Each loop is tied 
into a 3i knot and the Monte Carlo dynamics is such that each knot cannot translocate 
from its loop to the other one (Fig. 3). The two loops also remain unlinked. Since 
the number of configurations for the whole ring is maximum when one of the loops 
is much longer than the other one, most configurations break the symmetry between 
the two loops showing a marked length unbalance. Typically, in one of the two loops 
the knot has a very large share of the whole ring at its disposal, while the other loop 
is just long enough to host its knot. Consenquently, we choose to always identify the 
length of the shorter loop with the length, I, of a trefoil knot inserted in a SAP. Such 
/ is always sampled now in situations in which the two ends of the knot come close to 
each other in the neighborhood of the slip link. However, one can hope this statistic 
to be representative of more general samples if laws like Eq.(l) are considered. Indeed, 
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Figure 3. Sketch of a 3i knot (right) forced to its typical length by the entropic 
competition between knotted loops. To simulate the slip link we impose the constraint 
that two fixed parallel edges of the polygon have always to be kept one lattice distance 
apart 



for the average of this / the power law behaviour reported in Eq.(l) still holds, with 
t = 0.74 ± 0.05, consistent with the previous estimate (Fig. 4). In the case of two 
unknotted loops we verified that the average length of the shorter one does not grow 
appreciably with N, further supporting our interpretation of I as the knot length in the 
3i vs 3i case. 

This last method of length determination needs a suitable competition between the 
two loops [3U]. If only one of them were left unknotted, a too strong dominance of 
the knotted loop would result, making it almost impossible to sample configurations in 
which the longest loop is the unknotted one. Indeed, since the knot is rather localized, 
its freedom to place itself along the loop multiplies by a factor roughly proportional to 
the length of the loop itself the number of configurations accessible in its absence. By 
putting a 4i knot against another 4i knot, within error bars, we found the same value 
of t (Fig. 4). Thus, like in the previous strategy, t appears universal with respect to 
the prime knot type. We also carried out simulations in which one loop contained a 
trefoil, while the other one was tied into a composite knot made, e.g., by the product 
3i#3i of two trefoils In this case the loop hosting the composite knot almost 

constantly dominates in length because both components are rather localized and free 
to move independently along it. In spite of the different conditions of entropic balance, 
we verified that also in this case the average length of the (shorter) loop with a single 
knot satisfies Eq.(l) with the same t ~ 0.75 (Fig. 4). 

The above results have implications for the leading scaling correction to the singular 
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Figure 4. Log-log plot of the average value of the length I of the shortest loop as 
a function of N. Different curves correspond to the presence of different knots in the 
two loops. Best fit estimates of t are 0.74 ± 0.05, 0.75 ± 0.03, 0.77 ± 0.10, 0.74 ± 0.05 
for 3ivs3i, 3ivs3i#3i, 4ivs4i and 3ivs3i#3i#3i, respectively. 

behavior of quantities like the average radius of gyration of the SAP, (f? 2 ) 1 / 2 ~ 
AN V (1 + BN~ A ). Field theoretical renormalization group (RG) methods [31] and most 
Monte Carlo techniques allow to determine exponents like v ~ 0.588 and A ~ 0.50 
in ensembles in which configurations with all possible ring topologies are sampled, v 
is believed to remain the same in ensembles like those considered here, with fixed knot 
topology With our first cutting and closing method, we could check that the ring 
portion occupied by the knot has an average radius of gyration ~ l u °, with uq ~ 0.60 ~ v. 
This means that a length diverging as (Rq) ~ N vt exists in addition to the leading 
one (i? 2 ,) 1 / 2 ~ N u . Standard scaling arguments then lead to expect a possible strong 
correction exponent A = 1 — t ~ 0.25 for the ensemble with fixed knot topology. Such 
stronger correction could be quite difficult to identify numerically at relatively short 

7v mm. 

Our strategy based on entropic competition has the advantage of being applicable 
also to situations in which the polymer is compact. Indeed, since this method does 
not involve cutting and closing procedures, there is no risk that the knot topology gets 
altered, a highly probable event for compact configurations. This enabled us to address 
another key open issue in the field, namely the possible occurrence of knot derealization 
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in the compact regime ^H] • We included nearest neighbor attractive interactions for the 
two loop model and simulated it extensively at T ~ 0.53Tg, where Tg is the theta 
temperature (2HES1- Below the theta point it becomes more difficult to sample long 
SAP configurations. Indeed, at T < Tg, as the critical edge fugacity is approached 
from below, the grand-canonical average number of SAP edges, (N) g _ c _, undergoes a 
first order infinite jump from a finite value, rather than growing continuously to infinity 
as in the T > T@ case. In spite of this difficulty, we could get rather clear evidence 
that the canonical average (I) grows linearly with N up to iV ~ 700, implying t ~ 1, 
for the case of two competing 3i knots (Fig. 5). Like when dealing with the good 
solvent regime, in order to interpret the length of the shorter loop as the knot length, 
it is important to check what happens to the length of a loop without knot. We got 
evidence that in the compact regime an unknotted loop is also delocalized, i.e. shows 
an average length proportional to N. However, the proportionality constant is about 
one order of magnitude smaller than that observed for the knotted loop in Fig. 5. 
Thus, knot derealization seems to superimpose itself to a phenomenon of unknotted 
loop derealization which is qualitatively similar, although quantitatively much less 
pronounced. The entropic competition between loops in the compact phase does not 
lead to a strong prominence of one of them on the other, like at high T. To the contrary, 
all loop length ratios are realized with comparable probability. In this new scenario also 
the entropic mechanisms determining knot localization in the good solvent case cease 
to be active. 

The derealization of prime knots in compact polymer rings opens intriguing 
perspectives on how topological entanglement combines with more geometrical measures 
of complexity in polymers. One such measure for a fluctuating SAP is the mean 
number of crossings, N c , one obtains by projecting on all possible planes a given 
configuration [T^j. This quantity can be studied also for ideal knots ^H], and has been 
shown to be correlated with the mobility under electrophoresis of non-ideal knotted 
polymer rings with the same topology [3^. For ideal knots, N c is expected to grow as 
l^ 3 , where lu is the length of the knot (four third power law)[TT| H2]. For fluctuating, 
nonideal polymers one can consider the average of N c over all configurations. It is 
expected that quite generally, for a compact SAP, (N c ) ~ N 4 ^ 3 jHHj- This result should 
indeed hold even for open walks and be independend of topology for knotted SAP. The 
fact that we find (I) ~ N for a fluctuating, knotted compact SAP suggests that for real 
knotted polymer rings a statistical generalization of the four third power law of ideal 
knots should hold: (N c ) ~ (0 4//3 - I n the case of ideal knots the growth of 1^ corresponds 
to an increase of the topological complexity of the knot considered [ill El E3j • For real 
knotted rings the topology remains fixed with increasing N, while fluctuations are able 
to produce the same type of growth of (I) versus {N c ). 

In conclusion, we gave a consistent and robust statistical description of the 
localization properties of prime knots in polymers. Weak localization and derealization 
hold in the good solvent and collapsed regimes, respectively. 

This work was supported by MIUR-COFIN03, FIRB01 and INFM-PAIS02. 
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Figure 5. Plots of (I) of the shortest loop as a function of N. Both curves correspond 
to the presence of one 3i in each loop. The bottom curve refers to T >> Tq whereas 
the top to T < T e . 
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